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Abstract 

The well known relation for ideal classical gas Ae^ = kT^Cv which does not remain valid for quantum system 
is revisited. A new connection is established between energy fluctuation and specific heat for quantum gases, valid 
in the classical limit and the degenerate quantum regime as well. Most importantly the proposed Biswas-Mitra- 
Bhattacharyya (BMB) conjecture (Biswas et. al., J. Stat. Mech. P03013, 2015.) [T] relating hump in energy fluctuation 
and discontinuity of specific heat is proved and precised in this manuscript. 


1 Introduction 


An increasing attraction towards the subject of quantum gases is observed, after it was possible to create BEC in mag¬ 
netically trapped alkali gases [H 121 S] as well as experimental confirmation of Fermi degeneracy l^. Different theoretical 
and experimental studies analysing the effects of temperature dependence of energy and specific heat of ultracold Fermi 
gases[5l[6l[7], momentum distribution for harmonically trapped quantum gas[8l[9], temperature dependence of the chemi¬ 
cal potential[Tn], critical number of particles for the collapse of attractively interacting Bose gas[TT], Casimir effect [1^113] . 
equivalence of Bose and Fermi system[T31 [12], q deformed systemslUl [T7] have already been reported. Although a lot of 
theoretical studiesm m HOI HD HI m [Ml US are done on quantum gases trapped under generic power law potential, 
none of these contained detailed discussion on energy fluctuation of trapped quantum gases, until the recent paper of 
Biswas et. al. [1]. In the case of ideal classical gas, specific heat Cy is regarded as energy fluctuation Ae^ as the Ae^ is 
related to Cy as, Ae^ = kT'^Cy. But this relation becomes invalid for both types of quantum gases (Bose and Fermi) in 
the quantum degenerate regime for free quantum gases[T]. Biswas et. aL[T] have conjectured (BMB conjecture) a relation 
between the discontinuity of Cy and energy fluctuation. According to the BMB conjecture, the appearance of a hump in 
over its classical limit may indicate a discontinuity of Cy . They have shown this to be true for free and harmonically 
trapped Bose gases [I]. They have also mentioned without proving that the the inverse of this statement may not always 
be true. But this conjecture is yet to be proven for any quantum system trapped under generic power law potential in 
arbitrary dimension and is an open problem. In this manuscript, we have proved and precised BMB conjecture for ideal 
quantum gases trapped under generic power law potential, U = in. d dimension. Thus, in principle one 

can reconstruct the results Shyamal et. al., choosing all rii = 2 for harmonically trapped system and all rii = oo for free 
systems. Beside this, a relation is established between Cy and Ae^ which is valid not only in the classical limit but in 
the quantum limit as well. 


The manuscript is organized in the following way. In section 2, we have determined the grand potential in an uni¬ 
fied way for both types of quantum gases, from which we are able to calculate the quantities such as Cy and Ae^. In 
the next section we elaborate two theorems which eventually guide us to prove the conjecture. Results are discussed in 
section 4 and the paper is concluded in section 5. 

2 Grand potential, specific heat and energy fluctuations 


For a quantum gas, the average number of particles occupying the i-th eigenstate and the grand potential , are given by 
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where, a = 1(—1) stands for Bose system (Fermi system) and z is fugacity. Let us consider an ideal quantum system 
trapped in a generic power law potential in d dimensional space with a single particle Hamiltonian of the form, 

d 

€{p,Xi) = bp’- — (3) 

where, b, I, at, Ci, rii are all positive constants, p is the momentum and Xi is the i th component of coordinate of a particle. 
Here, Ci, Oi, rii determines the depth and confinement power of the potential and I being the kinematic parameter and 
Xi < Qi- As ll^-l < 1, the potential term goes to zero as all rii —oo. Using I = 2, b = -^ one can get the energy 
spectrum of the hamiltonian used in the literatures [nnniiiiiiii]- If one uses I = 1 and b = c one finds the hamiltonian 
of massless systems such as photons pH] . 

The density of states for such system is miis]. 


p(e)=C(6,U,')e^-i (4) 

where, C{b, Uj) is a constant depending on effective volume V);[I1]. For the detail derivation of density of states see Ref. 
[25| . Now, replacing the sum by integral we obtain the grand potential. 


Vi 

q = qo + j^Li^+i{cr) 

Note, Liq{m) is known as polylog function in the literature whose integral representation for Re{m) < 1 is[15) 
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It is a real valued function if m G K and — oo < to < 1. Also, the effective volume Uj, effective thermal wavelength 
along with y and cr are defined as. 
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z , Bose system 
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The detail idea of effective thermal wavelength and effective volume for trapped quantum gases can be found in[531 PH 
[25l [26] . But note, when I = 2 and b = ^ from Eq. (8) we get Aq = ^ which is the thermal wavelength of 

nonrelativistic massive fermions as well as massive bosons. However, it should be noted that, when I ^ 2, X' then depends 
on dimension. With d = 3 and d = 2, thermal wavelength of photons (boson) and neutrinos (fermion) are respectively]!^] 
2 t^i/ 2 kT ^’^'1 ( 27 r)i/ 2 fcr ^liich can be obtained from from Eq. (9) by choosing b = c, where c being the velocity of light. But 
one needs to consider the effects of antiparticles to calculate the thermodynamic quantities of ultrarelativistic quantum 
gas|27]. So, one can reproduce the thermal wavelength of both massive and massless particles from the definition of 
effective thermal wavelengtlQ X' with more general energy spectrum. It is also seen that the effective volume^ V/ is a 

^For detailed conceptual overview of effective thermal wavelength see Ref. m 

^Effective volume is also referred as pseudovolume, for detailed conceptual overview on how equation of state of trapped quantum system 
can be obtained see Ref. m 


2 










very salient feature of the trapped system as this allows us to treat the trapped quantum gases [TU [131 HH US] to be 
treated as a free one. Difference between and Vd is that, the former depends on temperature and power law exponent 
while the latter does not. But as all rn —s- oo, approaches Vd- The great benefit of evaluating DJ and A' is that 
they enable us to write all of the thermodynamic functions of the trapped quantum system to be expressed in a compact 
form similar to those of a free quantum gas [Ml [23]. It is well known that the Bose and Fermi functions do represent 
thermodynamics of Bose and Fermi system and can be written in terms of Polylogarithmic functions, 

= »(=) = Y(r)L 

Point to note, as z —>■ 1, the Bose function gx{z) approaches C(x): for X > llIH]- From the grand potential we can now 
determine, 


u. ATirj. ■^*x+i(o^ 

^ ^ = NkTx T - ! \ 
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In the classical limit of quantum gases, Cy equals to Nkx- Now, the energy fluctuation 

Ae^ = = X! = J dep{e)e'^n{e) - (J dep{e)en{e) f 


= (kTrixix + l) 
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Li\{a) 


(16) 


So, it is clear from Eq. (15) and (16) that, Ae^ ^ kT’^Cy- But it is valid in the classical limit as z —>■ 0, Eq. (15) and 
(16) deptics, 

Ae^ = kT'^Cv = Nx{kTf 


(17) 


But we can establish such a relation valid within the whole temperature range. Note, 
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where, in the last line we have used ^(f^)^= (§^)^ v' CS- temperature limit the second term of Eq. 

(18) becomes zero and Eq. (18) coincides with Eq. (17). It can be easily justified that equation (18) is valid not only in 
the classical limit but also in the quantum degenerate regime. 

Point to note that expression of Cy and Ae^ represents both Bose and Fermi system in a unified approach. In the 
case of Fermi system. 


Cy = Nk[x{x + 1)^%^ - X^T^] 

/xM) /x-iM) 

Ae^ = {kTf[x{x + 1)%# - X^%^ 


(19) 

( 20 ) 


The above equations coincides exactly with the results of Ref. [I] appropriate choice of Ui and d. Writing the expressions 
for Bose system (per particle). 


^[x(x + 1)f^5x+iW-X^ if^] ,T>re 

/cx(x + i)i)^C(x + i) .t<t, 

(fcT)2[x(x + l)«^-y2£|^] ^t>T, 

(fcT)2[yli(^)x-yl2(^)2x] " ,r<r. 
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Here Tq denotes critical temperature and, Ai and A 2 are defined as, 

,C(x + 2) 


^1 = x(x + l) 
^2 = [x 


C(x) 
C(X+ l)i2 


C(x) 


(24) 


Eq. (22) agrees with the expressions for free and harmonically trapped Bose system in three dimensional space[T]. It is 
noteworthy that, for both types of trapped quantum gases, Cy approaches its classical value x^K pS] , 


3 Theorems regarding Cy and Ae^ of Bose gas 

In this section we present the necessary theorems to prove the conjecture. As it was shown by Shyamal et. al. that 
a hump does exist in Ae^/kT'^Cy in the condensed phase for harmonically trapped Bose gas, we need to find a gen¬ 
eral criteria to locate a hump in Ae^/kT^Cy over its classical limit in arbitrary dimension with any trapping potential. 
As there is no condensed phase in ideal Fermi gas trapped under Dotential|29j. this theorem bears no significance for them. 


Theorem 4.1: Let an ideal Bose gas in an external potential, U = 
phase in Ae^/kT'^Cy over the classical limit if, 

Al 


i=i^ hump will exist in the eondensed 


X^2.3 


Proof: 

Re-writing Eq. (22) in the condensed phase, 


Ae2 

(fcT)2 

df 


= f^r) = Ait^ — A2T^^ 


/'(^) = ^ = ^iXT^-' - 2A2Xt^^-^ 


The condition for maximum, 
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The hump will be over its classical limit if, 


/(t) 


> X 

r=To 

X A J^X 


Ait^ - > X 

X< 1 (X) = A 


(27) 


From the Table 1, one can conclude relation (27) will be maintained (a hump will exist in the condensed phase over the 
classical limit) when, 

X > 2.3 (28) 


Theorem 4.2: Let an ideal Bose gas in an external potential, U = 
(i) Cy will be discontinuous atT = Tc if, 

d ^ 1 ^ 

— T -l- y — >2 

/ / ^ n • 
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Table 1: Listed values of 7 (x), generated using Mathematica (Correct upto four decimal points). 


X 

7(x) 

1.3 

0.8553 

1.4 

0.9756 

1.5 

1.1022 

1.6 

1.2350 

1.7 

1.3736 

1.8 

1.5181 

1.9 

1.6683 

2.0 

1.8240 

2.1 

1.9853 

2.2 

2.1519 

2.3 

2.3239 

2.4 

2.5010 

2.5 

2.6834 

2.6 

2.8709 


(ii) And the difference between the heat capacities at constant volume, a.t T = Tc as 

C(x) 


-Cv\. = Nkx^ 


C(x-1) 


Proof: 

As T ^ Tc, z ^ 1 and 7 —S' 0, where rj = — Inz. For T —)• T+, 

CviT+) =iVA:[x(x + l)^5x+iWL^i 

= Nkixix + l)l^C(x +1) - x^:r%L-.J (29) 

^ yx-iv^/ 

As the denominator of the second term of the right hand side contains g^-i{z), we can not simply substitute it with zeta 
function as z —>■ 1. So, using the representation of Bose function by Robinson 

( 30 ) 

' i^O 

we get from the above, 

CviT+) = Nkixix + l)^C(x + 1) - 

On the other hand 

Cv(Tc)=iVfc[x(x + l)^C(x + l) (32) 

Taking the difference between Cv{Tq) and Cv{Tq ), we get, 

ACy (33) 

Which dictates, ACy \t=Tc ^e non zero for x > 2 . So, Cy will be discontinuous when x > 2 and thus completing 
first part of the theorem. 
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As, X should be greater than two for ACy at T = Tc to be non-zero, we can re-write equation (21), by substitut¬ 
ing g^-i{z) by zeta function. 

Cv{T+) = Nk[x{x + 1)3^C(X + 1) - (34) 

From Eq. (32) and (34) we can write. 

ACy Cv -Cv (35) 

Note that, the same result is also derived by Yan et. al. for Bose gas trapped in symmetric power law notentiaipS]. 

Now, based on the above two theorems we can come to the conclusion. We have seen, a hump will exist in Ae^/fcT^Cy 
over the classical limit when y 2.3 and Cy will be discontinuous while y > 2. Therefore the existence of hump in 
Ae^/fcT^Cy over the classical limit automatically ensures discontinuity in Cy. But a discontinuity in Cy does not imply 
a hump in /kT^Cy over the classical limit when the value of y is in between, 2 < y < 2.3 and thus, proving the 
conjecture. 


4 Results and Discussion 

In this section we discuss energy fluctuation and specific heat in detail and check the prediction of the above theorems. We 
have illustrated in this paper how specific heat can differ from the energy fluctuation for the entire range of temperature. 
It is important to note that the difference between the probabilities of the classical and the quantum gas arises essentially 
from the nonzero fugacity of the quantum gas. 



T 

Figure I: Energy fluctuation ideal trapped Bose gas as a function of r = with different power law potentials. 


In case of trapped quantum gases all thermodynamic quantities are expressed by polylogarithmic functions depend¬ 
ing on fugacity and y. Thus apart from fugacity, the value of y bears the signature of difference between different 
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Table 2: Status of energy fluctuation and specific heat of Bose system trapped under generic power law potential. 


Range of y 

Hump over classical limit in Ae^ /kT'^Cy 

Discontinuity of Cy 

0 < y ^ 2 

2 < y < 2.3 
y ^ 2.3 

no hump over classical limit 
no hump over classical limit 
hump over classical limit 

continuous 

discontinuous 

discontinuous 


quantum systems. And as seen from the theorems the value of x dictates whether there will be a hump over the classical 
limit as well as the discontinuity of Cy- In figure 1 we have described the influence of different power law potentials on 
energy fluctuation of Bose system. It is clearly seen that, the Ae^/fcT^Cy has a hump way over its classical limit when 
X > 2.3. At X = 2.3. the hump is just over its classical limit. There is no hump over the classical limit when y < 2.3, All 
of these are in accordance with theorem 4.1. It is also noticed that, results in Shamyal et. al. [I] are also in agreement 
with the theorem. In their manuscript they found a hump over the classical limit in three dimensional harmonically 
trapped Bose system where x=| + | = 3> 2.3. Although they did find a hump in two dimensional harmonically 
trapped Bose system but this hump was below the classical limit. In this case, x=| + |= 2< 2.3 i.e., no hump over 
the classical limit. Therefore, it can be said that the theorem 4.1 can perfectly determine whether the humps will be 
below or above the classical limit. 

Now figure 2 illustrates Cy of Bose system with different trapping potentials. It is seen from the figure that, Cy 



X 


Figure 2: Specific heat of ideal trapped Bose gas as a function of r = with different power law potentials. 

is continuous when y ^ 2 and it becomes discontinuous when y > 2, in agreement with theorem 4.2. Now, as y ^ 2.3 
denotes a hump in Ae'^/kT^Cy over its classical limit, this automatically depicts discontinuity in Cy. Thus, we can 
conclude that the appearance of a hump in Ae^/fcT^Cy over its classical limit does indicate a discontinuity in Cy but a 
discontinuity in Cy does not conclude the appearance of a hump in Ae^/kT^Cy over its classical limit because disconti¬ 
nuity in Cy may arise even if 2 < y < 2.3 but no hump in Ae^ will exist in this interval of y. On the other hand, y ^ 2.3 
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will denote a discontinuity in Cy as well as the appearance of a hump in Ae^/fcT^Cy over its classical limit (see table 2). 


5 Conclusion 

In this manuscript we have restricted our study in the case of ideal quantum gases trapped under generic power law 
potential and proved the BMB conjecture for these types of systems. Point to note, as no hump in Ae^ or no discontinuity 
in Cv is noticed in ideal Fermi gases for any trapping potential. So, the theorems and the concluding relation between 
energy fluctuation and Cy remain significant for ideal Bose systems only. It will be interesting so see the status of the 
above theorem for interacting quantum systems. Also it will be very intriguing to generalize the theorems for relativistic 
quantum gases. 
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